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1 ABSTRACT 

\ ! 

This  report  describes  a method  for  obtaining  fine  frequency  resolu-  | 

tion.  The  technique  employs  the  partitioning  of  a large  data  sequence  | 

by  performing  small  size  FFTs  and  other  digital  signal  processing  tech-  ! 

niques  to  achieve  the  finer  resolution.  The  method  yields  approximate  j 

results.  The  accuracy  appears  to  be  dependent  on  Vernier  bandwidth,  | 

signal-to-noise  ratio  and  signal  spectrum.  A FORTRAN  program  is  available  i 

in  the  Appendix. 

H 
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GLOSSARY 


(n)j 

X (k) 

A f 

fs 

N 

M 

L 

FFT 


Time  data  sequence 

Discrete  Fourier  coefficient  at  K—  frequency 
Resolution  bandvddth  (in  Hertz) 

Sampling  frequency 
Time  sequence  length 
Number  of  partitions 
Vernier  bandwidth  (in  Hertz) 

Fast  Fourier  Transform 
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2.0  INTRODUCTION 

The  concept  of  the  ZOOM  FFT  was  first  suggested  to  the  authors  by 
Mr.  Stephen  Phanos,  Sperry  Gyroscope,  Great  Neck,  New  York,  fiach  of 
the  original  analysis  at  the  New  London  Laboratory  was  done  by 
Dr.  Nasir  Ahmed,  Kansas  State  University,  a summer  employee  at  rWSC, 

New  London. 

The  Zoom  technique  increases  the  sample  size  of  a vernier  band, 
by  creating  a new  narrowband  time  sequence  from  information  taken  from 
successive  large  size  but  course  resolution  FFTs.  Thus  finer  frequency 
resolution  is  obtained  that  is  an  approximation  of  the  true  spectral 
estimate  arrived  at  by  a Discrete  Fourier  transform  (DFT)  computation 
of  the  original  wide  band  time  sequence. 


MATHEMATICAL  DESCRIPTION 
Define  an  N point  real  sequence  (n^  , where 

= Xm,  X u"),  . ..  X<N-i) 


(3.1) 


The  DFT  (1)  yields, 


Tlie  resultant  frequency  resolutionAf  is  then  given  by: 


Af  = is 

N 


v/here  fs  is  the  sampling  frequency. 


(3.3) 


Consider  now  a new  sequence [>>j,  which  is  4 times  as  long  as  the  original 
sequence,  and  is  defined  by 


where  M = 4N. 


y(o),  •X(K-i') 
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The  choice  of  4 par Ci Cions  ls«  of  course*  arbitrary  and  used  for 
Illustrative  purposes.  In  the  general  case,  M is  equal  to  BN,  or 

M-RN 

where  R is  an  interger. 


By  taking  DFT's  of  the  4 partitions  each  of  length  N successively,  the 
resultant  four  sequences  are  defined  as; 

[Xiw}=  

* XiC°\  XjCi) Xz  w-i") 

{XjMj  = X3  W,  X-5 

PX/(  i . X<(  (3), . . - -X-1  c^-  ^ ) 

The  rationale  for  partitioning  is  twofold.  First,  the  available 
FFT  hardware  may  be  constrained  to  do  only  an  N point  DFT.  Second,  in 
order  to  perform  the  computations  In  real  time  it  may  be  desirable  to 
begin  processing  each  partition  as  it  becomes  available  or  in  parallel 
with  other  operations  where  each  sequence  has  a frequency  resolution 
equal  to  A ^ • 


tk>w  define  a limited  frequency  band  consisting  of  L frequency 
points,  and  using  the  same  band  for  each  of  the  4 previously  defined 
sequences,  write  these  sequences  as; 


voKe.vc.  K=P,F+-1- 
else  uiVere 


p-^L- J- 

(2,7) 
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The  inverse  FFT  of  each  successive  L point  sequence  is  computed. 
Even  and  odd  symnetry  is  forced  on  the  real  and  imaginary  Fourier 
coefficients  such  that  the  resultant  time  sequence  is  real.  The 
following  relationship  is  used: 

P-H,-!- 


Y_  X^U')EXPrJ_^"~ni. 

JL‘? 


and  r = 1, 2,3,4 


v\=  o.a. . - - - 


(3.8) 


Then  we  juxtapose  the  resultant  time  sequences.  Defining  this  new 
sequence  as 


CJ.'t) 


and  by  taking  the  OFT  of  this  new  sequence  according  to 

H-L-l. 

Xc-V-  C310) 


V^=o,i,  — 


results  in  a frequency  resolution  given  by 


A f = fs 
4M 


(311) 


or  4 times  as  fine  as  the  original  resolution.  In  effect,  the  Zoom  FFT 
has  performed  a large  size  FFT,  for  a Vernier  bandwidth  v;ith  resolution  4f ! 

A more  rigourous  mathematical  presentation  can  now  be  offered  to 
illustrate  the  frequency  approximation  process. 

Define  an  N point  time  sequence,  and  partition  the  sequence  into  M 
partitions,  the  resultant  M FFT' s can  be  represented  by 
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N\ 


XC>n-+KM')  EXP  "IHLilii 


where  maO,  1,  ...,  M-  1 
n a 0,  1,  N - 1 


iUl] 


which  yields  M frequency  sequences  of  size  N - 1.  Discarding  certain 

M 

frequency  components  by  selecting  X (R)  where  R is  defined  as 

rr*\ 


V\  ^ V\  u>V\e.vc- 

*^1.—  ~X\^=a  L 

R=  o,  1.,  - - - . (.u-JL') 


(113) 


This  relation  yields  M frequency  sequences  each  of  size  L.  Taking  the 
inverse  transform  of  each  X (R)  by 


U-X 

Rsn  b=o, 


(31i)) 


a.  - - - (.L-i'i 


which  yields  a M point  time  sequences  each  of  size  L.  By  forming  a new 
time  sequence  of  size  ML  by  juxtaposing  in  time  according  to 


- a. 

X(j)=  ^y^cb) 


where  d = b + ML  for 

b a 0.  1 - - (L-1); 
maO,  1----  (n-1) 
and  d a 0,  1 . . . (L-1)  . . (ML) 

and  taking  the  FFT  of  the  time  sequence  of  size  ML 


(315; 

(3Hj 


Xw-  ^W«e*pf-1^]  (317) 


S = 0,  1 . . . (ML-1) 
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Equation  (3.17)  yields  a frequency  sequence  of  size  ML, 


Given  equations  3,7  through  3.8,  the  Zoom  FFT  process  can  be  written  as: 

(3.16) 


/»t«o 


By  setting  the  E equal  to  N-1,  the  degenerate  case,  it  can  be  seen  that 


M 


equation  (3.18)  yields  the  following  equation  (2)  for  a partitioned  OFT 
case  namely 


-0, 


M-1  ^ 

^EXP(-32jI2^x')  ^XP 

IC^D 


/'r=o 


(5. 1'fj 

Thus  the  only  time  the  Zoom  FFT  is  equivalent  to  the  DFT  is  for  the 
degenerate  case.  ^Viy  lesser  size  of  E,  other  than  N-1,  yields  an 
eipproxiraate  answer.  M 


In  the  limit,  as  — ^^Ii“l)>  Che  approximation  goes  to  the  exact 

M 

answer.  The  distortion  in  the  algorithm  is  introduced  by  attempting 
perfect  filtering  in  the  frequency  domain,  thus  by  increasing  the  filter 
bandwidth,  less  error  is  introduced  in  the  time  domain. 


Examples  of  the  Zoom  FFT 

To  present  the  effect  of  changing  the  inverse  FFT  v^indow  (L),  two 
synthetic  tones  v;ere  digitally  produced,  one  at  430  Hz  and  another  at 
401  Hz,  the  latter  10  dB  lov;er  than  the  former.  A DFT  and  Zoom  FFT 
was  performed  on  this  date. 


Figure  3.1  presents  the  plot  of  the  DFT  case  using  an  FFT  algorithm 
(2’’,  with  1 Hz  resolution.  This  output  was  used  as  a reference  for 
comparison  purposes.  The  Zoom  FFT  was  tnen  performed  on  the  same  time 
data  sequence  of  1024  points. 
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First  the  sequence  was  partitioned  into  four  smaller  data  sequence 
of  256  points  each.  Next  an  FFT  v;as  taken  of  each  sequence;  this  FFT 
had  a resolution  of  4 Hz.  This  frequency  domain  data  was  then  used  as 
an  input  to  the  new  Zoom  algorithm  to  obtain  1 Hz  resolution. 

The  vernier  band  of  interest  was  selected  to  provide  a vernier  for 
the  401  Hz  tone.  Vernier  bandwidths  of  32  Hz,  64  Hz,  128  Hz,  and  256  Hz 
were  examined  and  shov/n  in  Figures  3.2  through  3.5.  Tlie  degree  of 
^proximation  can  be  compared  with  the  DFT,  quantitatively,  in  Appendix  B. 

Upon  examining  these  values  and  figures,  it  is  clear  that  as  the 
Vernier  window  is  increased,  the  error  of  the  approximation  to  the  DIT 
is  less. 

To  prove  this  conclusively,  the  window  was  made  equivalent  to  the 
partition  DFT  and  the  Zoom  FFT  was  performed.  The  results  are  shown  on 
Figure  3.6  and  are  exactly  equivalent  to  the  DFT  reference  case.  The 
results  substantiate  the  mathematical  interpertation  derived  in  the 
previous  section. 

An  interesting  aspect  of  the  attempt  a perfect  filtering  in  the 
frequency  domain,  is  the  possibility  of  reducing  the  resultant  side 
lobe  structure  by  Cosine  windowing  (4).  This  window  was  applied  to  the 
real  and  Imaginary  frequency  sequences  before  the  inverse  FFT  was  taken. 
The  degenerate  case  (v/here  the  ^om  FFT  is  equivalent  to  the  partitioned 
DFT)  was  examined.  Tlie  results  for  one,  two,  and  four  partitions  are 
shown  in  Figures  3.7  through  3.9.  It  is  apparent  that  windowing  in  the 
frequency  domain  in  this  manner,  increases  the  amplitude  of  the  side 
lobe  structure;  the  amplitude  also  appears  to  increase  with  the  number 
of  partitions  used  in  operating  the  original  time  sequence. 

Figures  3.10  through  3.12  shows  DFT  and  the  Zoom  FFT,  with  vernier 
bandwidths  of  64  Hz  and  128  Hz  using  the  two  tone  case  plus  additive 
random  noise  with  a mean  10  dB  below  the  401  Hz  tone.  As  expected  the 
Zoom  FFT  performed  best  for  the  larger  vernier  bandwidth. 


4.0 


SUICtARY 


The  Zoom  FFT  algorithm  can  be  used  as  a Vernier  FFT  algorithm  with 
approximate  results.  This  approximation  approaches  to  a partition  DFT 
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bandwidth  increases.  Cosine  windowing  in  the  fre- 
to  yield  negative  results. 


CONCLUSION 


Tlie  Zoom  FFT  appears  to  be  applicable  in  cases  where  a Vernier  FFT 
algorithm  is  needed  and  a real  time  constraint  is  such  that  other  (2.5) 
Vernier  algorithms  cannot  be  Implemented.  By  degrading  the  results,  the 
Zoom  FFT  can  be  implemented  using  very  small  FFT's,  thus  saving  appreciable 
time.  If  the  time  constraint  is  such  that  the  desired  frequency  resolution 
can  only  be  attained  by  reducing  the  Vernier  bandwidth,  the  accuracy  will 
degenerate  rapidly  in  cases  of  low  signal-to-noise  ratio  and  closely  spaced 
spectral  lines. 
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C ZoONi  FF'T  **■*.  j.KtHiilE  , C.NAWKOCKI  ♦**  JANUARY  3»1973, 


• OPlHATING  liJifAjCl  ions  ZOOM  FFT 


GaaC  1 MUot  KLPt-.AT£.D  FOR  EACH  DATA  SAMPLE 


L 


C C«kD 

lOLu.-.nS 

F OR  <iAT 

Ak&u'-'EN  T 

C i 

1-10 

llO 

,4NN 

- 

NUM6ER  OF  input  VALUES 

TO  BE  processed 

c 

ll-«iC 

ao 

ISR 

- 

integer  Sampling  rate 

c 

21-jO 

ao 

NpAR 

- 

number  of  partitions 

L 

31-*f0 

ay 

LFoiri 

- 

LO^.ER  FREOUEF'ICY  BIN  OF 

INTEREST 

L. 

41-jO 

ao 

UFdIN 

- 

HIGHER  FREiiiUENCY  BIN  OF 

INTEREST 

L 

5l-o0 

t 10. b 

F3«SE 

- 

fkeouency  of  vAfuable  tone  in  HE'RTZ 

L 

61  -f)0 

tiLA  jK 

w 

L L A'  b T 

1-..  fl 

^LA•A 

C ♦ ♦ * • 

^ « 4 

♦ ♦ 

+ t * ♦ 

3$. 

♦ j*;******#’*** 

* * * m * * * 

L 

u 

L SPEliFIlaT  lO.j  A jO  type  br"TEMENr5 


uiMENSiON  AX(lLi<i'+)  ,>Y(l(i24)  ,ZX(1024)  »ZY(1024) 

L)lME(4biON  PHl  (1024)  ‘ 

uIMENSiON  1 iTLt  a(3)  ,T  ITl.EY(5) 

LOMMuN  AMuuE^>lt00)  I 

uIwEnSiOI-^  fii''iE(ti) 
liNfEoEtA  UFijI'J 

tEFlriL  F-  (A  i tXK)  = II  ♦ivK  ♦ ( ( l-blGUl  I • 11*K,; ) )/<  j ♦ 54359738367 

C I 

^ i 

L 1MIUAL1..U  P<OG^Av  U.wITb  ANJ  IGb  ^'LOirP:R  i 

c 

llAtAi-AJ^i 

JF’kNTR:*! 

LAlL  MuUE  ( ANio  jES  , ij  ) 

KaSL  = i 

c 

c 

C jET  oF’  Ou  hA.'oE  A'iD  a*Y  IITl-ES  op  i‘i_uTj 

L 

IjD  -lAXTiO  , 0 
liiLtA(l)-'  Frt-0* 

1 1 |LEX(2)-'  )cNLY  ' 

HILL  X(3)-'  (H^.)  ' 

U 1L(  r n ) = 'Dp.//l  ’ 

U iLf.  r (2)  A WwL  1 >• 

T1  iLi.  Y(3)-’0./m2  ' 


.-1 
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ioO  READ(InCAK0»102)  NimN,  ISR,NPAR»LFBIn»UFBIN,F8ASE 
102  FORMAT(5I10fFlO,5) 

lFlNNN.6T,i024)  STOP  NNN 
IF(NnN.EQ.O)  go  to  900 

»*RlTE(iPRNrR»104)NNNi  iSR»NPAR»LFBlN,UFBIN,FBASE 
104  F0RMAT(1H1.///ICX»  *NNN  = » » If>r  5X , ’ I SR  = ' . 18 . 5X , * NPAR  =» . I5,5X,  ♦LFBI 
IN  =’.Ib,5X,«UFBlN  =' » I5,5X, *FBASE  = » , FI  0 . 5 » /// ) 

calculate  constants 

DF=FLOAT(lSR)/rLOAi iMNN) 

D1=1.0/FLuAT(ISk) 

DFF=DF*FLOaT(NPAK) 

COnST=DT/PuOAT (NnN) 

1SpAr=nNN/nPAR 
lDlFF=uFBil(“LFc3iN 
1DIFFi  = ID1fF-*1 
1D1FF2=2*IuIFF 
NRESilOIFFd*NPAK 
NDlS=Nr(E5/c  + l 
FuU^=Fi_OA|  (CFB1n-1  ) ♦DFF 
FHlGH=FLOA  + NUlb*DF-l  .0 
FIuIFF=1,0/FlOAT (1j1FF2) 

aR1TL( IPRNTRtllO)  iSPaR, 1D1FF2.NR£S*N0IS»FL0W»FHIGH,DF 
liO  FOKMAT(/»iOX» ’ ISPAr  =’ , I5»4X» • IOIFF2  » Ib»4X, *NRES  =* » I5,4X» 'NDIS 
i =.*  » I5,4X. 'FLUrt  =»  ,F10.3,4x, 'FHIGH  = • ,F10.  J,4X»  »DF  =•  ,F12.8.///) 
IF(lDlFF2.UT.lbPAR)  STOP  DIF2 


GENckATE  iNPuT  data 

*KITE(  IPRNTR»l.iO) 

120  format (///lox, ' Input  data'.///) 

Pi-3.l4l59c.6b4 

OAiNzSoRT liO.O) 

F Ki.Oc=430 . 0 

^^=5♦  ♦ I'j 

Ii=5?bl 

0MlN2=bQRl l 10,0) 

uU  IbV  jri.NNN 

aX  (U)=bIN(.c.O*Pl*F;jA5E*OT*FLOAT(  J)  ) 
TLMPrbAlN*bIi<i<-.0*Pi*FKLO2*FLOAT  ( J)  ♦DD 
aX(J)=xX( J)+TEmP 
li=F(U.KX) 

\/  = T1NO«M(FuOAT  {lI)/o4  3b973b36/.  ,S126) 

GO  To  126 
l*i6  COnTINjE 


AKlTE(iPR..|R.lt/) 

127  FORMaT(1H  , 'I’RotiLE  l*  ) 
l26  aX(J)=xX(  j)+V*gaIN^; 

y r(U)=o.o 

<1X(J)=aX(o) 

ZY(J)=YY (D) 
lb7  oOnTINuE 
u1uPRT=200 

iF  lN[.N.LT,i.TOpt-((  ) nTOPR  r=f'<i/l 
*RiTc( IPRNTR.2^4)  (XX(J) .YY(J) . J . J=1 . NTOPR T) 
,,c.4  F orMA  T ( b ( .if  9 . 5 . i 4 . 4X  ) ) 

A- 2 
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COSINE  WINDOW  INPUT  DATA 

CAUU  C05MTH(XX»NNN) 

IF|ICASC.*Ta)_60  10230  > 

C AlI' COSH^TH  rzV'r  NNN ) 


COMPUTE  UFT  OF  SIGNAL  TO  CHECK  ZOOM  FFT  ALGORITHM 

calL  first 

call  FFT(ZX»ZY»NNN,NNN,NNN»-I) 
call  SLCONU(TIME) 
wHlTE(IPRi'tTR»225) 

22b  F0rMaT(////, * OFT  check  CASE’) 
wRITE(IPRNTR»242) 

ND2PI=NNN/2+I 
DO  226  j=i,ND2Pi 

Pnl (J)=C0NsT*iZX<J)**2+ZY(J)**2) 

FREQ=Fl0AT(J-1)*DF 
uD=10.U*AuuGlO(MAX(PHl (j) ,i.0E-30) ) 

»KlTt(IPRNTR»246)  ZX(J) ,ZY( J) » J,FR£Q»pHl ( J) ,DB 
PHI(J)=0B 
226  cOlgTlNuE 

WKlTE.(iPrtNTR,344)  time 

CAlL  LINPlT (0.0.DBmIN,FREQ»DBMAX»PhI»1»ND2P1,TITLEX,TITLEY»0) 
call  PAGEo(AMOUES»Of 1»1) 

I4ASl=2 
230  continue 


CALCULATE  ZOOM  FFt 

call  First 

DO  240  K=X*NPAk 
1SpAr1=ISpaR*(k-1) 

1SpAr2=idIFF*;*(k-1) 
iNuExl=ISPARl4i 

lNuEX2=lSPARl+uFdIu-I 
C FORWARD  FFT  uF  INOIYiDUaL  PAkTITIONS 

call  FFT(aa(INUEX1) ,Yy(iNDtXl) ,ISPaR»ISPAR»ISPAR,-1) 

C PICK  OLI  desired  DINS  AND  TAKE  INVER:>E  FFT 

DU  234  L=i»IOlFH 
ZX(L)=xX(I(4DEx2+L) 

ZY  iL)zYYUwDEX2+L) 
fc34  COwTINjE 

C (jENLKATE  uVEr.  and  odd  SYMMETRY  BEFOKt  INVERSE  FFT 

ZY(I0IFF1)=0,0 
DO  23B  L=2»IDIFF 
NEO=IDlFF2-L+2 
ZX(NEG)=ZX(L) 

ZY  (NF.G)=-ZY(L) 
c-bb  uOf^TlNuE 

call  FFT(ZX»ZY» ID1FF2. IDIFF2»IU1FF2»+1) 

uO  239  I=1iIUIFF2 
aX(1sPaR24.I)=ZX(  I)  tFIDIFF 
I tYllSPAR2+l)=ZY(I)*FIDlFF 

' Zb9  continue 

' 240  continue  A-3 


o r--  o n r,  o o r>  rt  r 


C 


IM  NO.’ 

forward  FFT  of  LFBIN  to  UFBIN  bins  SA2201-584-72 

CAlU  FFT(XX»YY»NREs»NRES»NRES,-1) 

CALt  SECONO(TIME) 


display  kESUlTS 
wRITEUPRNTR*242) 

<h2  format (///icx» 'Complex  coefficient^' »//12x. 'real* *16x» *imag* »i3X» * 

IN’ »7Xf 'FREo*  »11X» 'POWER' ,10X» 'DB' »/) 

DO  250  J=1,NDIS 

PHi  (o)=COlMbT*  iXX(  J)**2+YY(  J)**2) 

FREQ=FlOW+FLOAT { j-i ) *DF 
Ob=l0.U*ALOS10(MAX(PHl (j) ,l.0E-30) ) 
write (1PRNIR» 240)  XX( J) ,YY( J) , J,FREO»PHI ( J) ,0B 
246  FOKMAT(2F2i),b,llOoFl5.8) 

PHI(J)=DB 
250  continue 

wKITEc IpRNTR.344)  TIME 

344  KirMaH/Z/.IOX, 'EXECUTION  TIMl  '»2nfa,'  SECONDS') 

CAuL  LINPi_f (PLOW,DjMlN,FHlbH»UBMAX,PHI,l,NuIS.TITLEX,TlTLEY»0) 

CmlL  PAGEb ( AVOutSf 0 » I » 1 ) 


oO  PACK  FOR  uEXT  CASE 
bU  TO  100 


lERi-iUjAjE  PRybRAM 

bUO  call  ExITbiAVOUES) 

STOP  200M 
ENu 


c ...  SUoROUUNi.  C.iMTH  ...  NOR-«ALIiED  TI«E  COSINE  -IN0O,  ..H.NUTT.LL 

SOoROUflNt  COSi“'.TtUxX»NNiM) 
u1i*iEnSiON  aX  1 1 ) 

T=0.?B3183jP7/FlOAT  iM.’iN) 

kETUKN 

cNu 
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TM  iNO. 

SA2201-584-72 


APPENDIX  B ' 
DFT  CASE 


391.00000000 

392.00000000 

393.00000000 

394.00000000 

395.00000000 

396.00000000 

397.00000000 


.00000000  -146,136652ao 
.00000000  -146.12798309 
,00000000  -145,08983803 
.00000000  -132,95911789 
,00000000  -131.25646591 
.00000000  -137,40895081 
.00000000  -123.81223202 


396.00000000 

399.00000000 

400.00000000 

401.00000000 

402.00000000 

403.00000000 

404.00000000 
405.0000COOO 

406.00000000 

407.00000000 

408.00000000 

409.00000000 

410.00000000 

411.00000000 

412.00000000 

413.00000000 

414.00000000 

415.00000000 

416.00000000 

417.00000000 

418.00000000 

419.00000000 

420.00000000 

421.00000000 

422.00000000 

423.00000000 

424.00000000 

425.00000000 

426.00000000 

427.00000000 

428.00000000 

429.00000000 

430.00000000 

431.00000000 

432.00000000 

433.00000000 

434.00000000 

435.00000000 

436.00000000 

437.00000000 

438.00000000 
43P. 00000000 
440.0000''000 


,00000000 
.00000000 
,04166749 
.16666679 
.04166590 
.00000000 
,00000000 
,00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
.00000000 
,00000000 
.00000000 
.00000000 
,00000000 
,00000000 
.00000000 
,00000000 
,00000000 
,00000000 
,00000000 
,00000000 
,00000000 
,00000000 
.00000000 
.00000000 
.00000000 
.41667493 
1 . 66666658 
.41665835 
,00000000 
.00000000 
,00000000 
,00000000 
.00000000 
,00000000 
,00000000 
,00000000 
.00000000 


-120,45800877 
-121.02896309 
-13,80202663 
-7,78150934 
-13,80219221 
-124,35157871 
-133,65653946 
-130.35192490 
-131,80841446 
-134,32049370 
-143.15020752 
-142,63742447 
-134,74664497 
-132,38457108 
-138,01848412 
-131.02366638 
-125,16299820 
-131,05418968 
-138,69378090 
-136,85265732 
-133.30390549 
-131,73151016 
-133,18720627 
-132.74911880 
-123,90300465 
-128,11253548 
-130,92435265 
-124,33635426 
-120,85968781 
-120,61209965 
-121.37391567 
-3,80202630 
2,21848726 
-3.80219907 
-121.91281891 
-120.18064499 
-121.98831367 
-12o. 57015610 
-129.03431892 
-135,20573616 
-124,93441391 
-130,26271057 
-135.79203606 


1 

1 

i 


I 

I 

i 
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TM  NO. 

SA2201-584-72 


ZOOM  FFT  WITH  VERNIER  BANDWIDTH 
OF  16  HZ;  1 HZ  RESOLUTION 


FREQ 

59P.00C00000 

393.00C00000 

394.00000000 

395.00000000 

396.00000000 

397.00000000 
39B.OOOOOOOO 

394.00000000 

400.00000000 

*♦01.00000000 

402.00000000 

403.00000000 

404.00000000 

405.00000000 

406.00000000 
4O7.0O00COOO 
40H.OOOOCOOO 


POWER 

.00000000 

.00060436 

.00066398 

.00038748 

.00000000 

.00087759 

.00027502 

.00010563 

.04166749 

.14367197 

.04155525 

.00005242 

.00000000 

.00281936 

.00221985 

,00084729 

.00000000 


DB 

-148.22646713 
-32.18701506 
-31,77844524 
-34.11746168 
-137,35529900 
-30.56709361 
-35.60634804 
-39.76223183 
-13,80202675 
-6,42627954 
-13.81374073 
-42.80472851 
-133.67460251 
-26.49850011 
-26.53676772 
-30.71970296 
-14b. 15300941 


m 


m 


TM  MO. 

SA2201-584-72 


ZOOM  FFT  WITH  VERNIZR  BANDWIDTH 
OF  32  HZ;  1 HZ  RESOLUTION 


FKtQ 


POWER  DB 


O84.0O00C000 

355.00000000 

386.00000000 
387«OoOOOOOO 

380.00000000 

389.00000000 

390.00000000 
39l,OoOQOOOO 

392.00000000 

393.00000000 

394.00000000 

395.00000000 

396.00000000 

397.00000000 

398.00000000 

399.00000000 

400.00000000 

401.00000000 

402.00000000 
403.0UOOOOOO 

404.00000000 

405.00000000 
406.0000C000 

407.00000000 

408.00000000 
40P,0u000000 
410.0UOOOOOO 

411.00000000 

412.00000000 

413.00000000 

414.00000000 

415.00000000 

416.00000000 


.00000000 

.00005125 

.00293665 

.00021421 

.00000000 

.00005102 

.00240252 

,00015874 

.00000000 

,00007590 

.00204946 

.00015091 

.00000000 

.00012189 

.00173565 

,00016308 

.04166749 

.15651207 

.05712260 

.00018715 

.00000000 

.00035623 

.00095036 

.00022076 

.00000000 

,00069765 

.00046813 

.00027215 

.00000000 

.00180507 

.00261760 

.00117454 

.00000000 


-143,58618164 
-42.90313005 
-25.32147932 
-36.69154072 
-153,53866386 
-42.92224503 
-26.19332552 
-37.99326468 
-146,06138992 
-41,19755411 
-26.88360810 
-38.21289301 
-137.40611649 
-39.14017534 
-27.60538626 
-37.87609577 
-13,80202663 
-6,05452156 
-12.43192029 
-37.27802038 
-133,69916153 
-34,48267174 
-30.22111678 
-36,56080675 
-143,19004631 
-31,56365204 
-33.29635000 
-35.65199280 
-137.96994781 
-27.43505621 
-25.62096124 
-29,30130863 
-14;?. 34171867 
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TM  NO. 

SA2201-584-72 


ZOOM  FFT  WITH  VERNIER  BANDWIDTH 
OF  64  HZ;  1 HZ  RESOLUTION 


375.00000000 

376.00000000 

377.00000000 

578.00000000 

379.00000000 

330.00000000 

381.00000000 

382.00000000 

363.00000000 

384.00000000 

385.00000000 

386.00000000 

387.00000000 

388.00000000 

389.00000000 

390.00000000 

391.00000000 

392.00000000 

393.00000000 

394.00000000 

395.00000000 

396.00000000 

397.00000000 

390.00000000 

399.00000000 
400.0000C000 
'+01.00000000 

402.00000000 

403.00000000 

404.00000000 

405.00000000 

406.00000000 

407.00000000 
40P.OOOOOOOO 

409.00000000 
410.0000COOO 

411.00000000 

412.00000000 

413.00000000 

414.00000000 

415.00000000 

416.00000000 

417.00000000 
418,00U0'?000 

419.00000000 


,00074553 

.00000000 

.00031923 

.00906725 

.00078177 

,00000000 

,00029906 

.00905400 

.00083644 

.00000000 

.00028903 

,00924582 

.00091046 

.00000000 

.00028675 

.00963630 

.00100819 

.00000000 

.00029138 

.01024772 

.00113724 

.00000000 

.00030311 

.01113142 

.00131020 

.04166749 

.17847740 

.09028197 

.00154734 

.00000000 

.00035306 

.0141J954 

.00186445 

.OOOOuOOO 

.00039742 

.01669311 

.00236694 

.OOOOuOOO 

.00046325 

.02055604 

.00316809 

.00000000 

,0005o520 

.02681831 

.00459815 


-31.27533031 
-144.56781960 
-34.95895863 
-20.42524529 
-31.06916693 
-146.43681908 
-35.24242973 
-20.43159628 
-30.77565813 
-143.71667480 
-35,39052916 
-20.34054399 
-30,40737224 
-154.66052818 
-35.42501068 
-20,16089869 
-29,96458721 
-146.61285400 
-35,35540724 
-19.89372516 
-29.44146633 
-137.52787018 
-35.18393993 
-19.53449607 
-26.82661724 
-13.80202675 
-7.40416781 
-10.07526124 
-26.10414982 
-133,78762436 
-34.52140199 
-18,49564648 
-27,24815679 
-143,04985237 
-34,00747442 
-17,77462769 
-2b, 22159362 
-136,04380608 
-33,34161929 
—lb, 87060618 
-24.96469426 
-130,77390671 
-32.47796774 
-16.71568513 
-23.37416792 
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TH  NO. 

SA2201-584-72 


ZOOM  FFT  WITH  VERNIER  BANDWIDTH 
OF  128  HZ;  1 HZ  RESOLUTION 


392.00000000 

393.00000000 

394.00000000 

395.00000000 

396.00000000 

397.00000000 

398.00000000 

399.00000000 

400.00000000 
40l.0(j000000 

402.00000000 

403.00000000 
404.0O0OC000 

405.00000000 

406.00000000 

407.00000000 
40B.00000000 

409.00000000 

410.00000000 

411.00000000 

412.00000000 

413.00000000 
hIV.OoOOOOOO 

415.00000000 

416.00000000 

417.00000000 

418.00000000 
419,no000000 

420.00000000 
421 .OuOOCOOO 

422.00000000 
425.0UOOOOOO 
424.OUO0OCOO 

425.00000000 

426.00000000 

427.00000000 

428.00000000 

429.00000000 

430.00000000 

431.00000000 

432.00000000 

433.00000000 

434.00000000 

435.00000000 
43ej. 00000000 

437.00000000 

438.00000000 

439.00000000 

440.00000000 


,00000000 

,00000633 

,00017162 

•00003836 

•UOOOOOOO 

.00000655 

,00017081 

.00003879 

.04166749 

.16495706 

.04677598 

.00003936 

.00000000 

.00000724 

.00017198 

.00004013 

.00000000 

.00000774 

.00017394 

.00004108 

.00000000 

.00000837 

.00017682 

.00004224 

.00000000 

.00000917 

.00018067 

.00004362 

.00000000 

.00001019 

.00016555 

.00004525 

.00000000 

.00001148 

.00019158 

.00004717 

.00000000 

.42011832 

1.63444418 

.40855611 

.00000000 

.00001527 

.00020800 

.00005199 

.00000000 

.00001609 

.00021912 

.00005502 

.00000000 


-146.16699600 

-51,98672581 

-37,65429115 

-44,16107035 

-137,40612030 

-51,83911848 

-37,67474937 

-44.11311065 

-13.00202607 

-7.82629091 
-13.29977119 
-44,04945660 
-133,66323853 
-51,40287304 
-37.64523792 
-43,96532536 
-143.28040886 
-51,11419535 
-37,59594727 
-43.86320400 
-137,98767090 
-50,77159023 
-37.52461386 
-43.74276209 
-139.13274765 
-50,37527418 
-37,43108559 
-43.60283470 
-133.16700281 
-49,91934681 
-37.31549263 
-43.44398117 
-130.86851501 
-49,40131044 
-37,17651129 
-43.26312876 
-121,38550854 
-3.76628375 
2.13370091 
-3.80746163 
-121,93789768 
-46,16065788 
-B6,81927Q19 
-42.04105206 
-129.02977180 
-47.42601252 
-3o. 59309101 
-4^,59447956 
-115.02561111 
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TM  NO. 

SA2201-584-72 


^\PPSNDU  C 

It  is  apparent  from  the  report,  that  for  a fixed  FFT  and  a desired 
resolution,  the  largest  Vernier  bandwidth  possible  should  be  used  to 
improve  accuracy.  If  the  algorithm  is  to  be  exercised  in  real  time, 
a constraing  to  the  bandwidth  is  made  by  the  hardware  multiply  times. 

-1 

We  can  model  the  FFT  and  FFT  respectively;  by  [C-lJ 
LMlnL  + LMlnLM  = Tm  (c-lj 

''Aiere  L is  the  vernier  bandwidth 

M is  the  number  of  juxaposed  sequences 
and  Tm  is  the  total  number  of  multiplies 

The  above  equation  is  unseperable  by  ordinary  techniques,  however 
we  allow  M to  take  on  interger  values  (M  = 1,  2,  4,  8,  16,  32)  and  increment 
L,  so  that  the  product  of  LM^024  (where  1024  is  a typical  hardware  FFT 
size),  constant  multiply  lines  can  be  drawn,  as  shown  in  Figure  C-1. 

Thus,  if  the  per  multiply  time  of  the  processor  was  a,  and  a 8 times 
finer  resolution  v;ere  desired,  by  multiplying  a times  point  picked  from 
the  constant  multiply  so  that  the  real  time  constraint  wa-;  ratisfactory 
and  the  largest  modulo  -2  L was  selected. 
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